We use a numerical method to obtain the weak coupling perturbative coefficients of local operators with lattice regularization. Such a method allows us to extend the perturbative expansions obtained so far by analytical Feynman diagrams calculations. In SU(3) lattice gauge theory in four dimensions we compute the first eight coefficients of the expectation value of the Wilson loop on the elementary plaquette which is related to the gluon condensate. The computed eight coefficients grow with the order much faster than predicted by the presence of the infrared renormalon associated to the dimension of the gluon condensate. However the renormalon behaviour for large order is quite well reproduced if one considers the expansion coefficients in a new coupling related to the lattice coupling by large perturbative corrections. This is expected since the lattice and continuum Lambda scales differ by almost two orders of magnitude.